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THE SOLUTIONS OF DIFFERENTIAL EQUATIONS OF THE FIRST 
ORDER AS FUNCTIONS OF THEIR INITIAL VALUES 

By G. A. Bliss 

The proof of the existence of solutions of the differential equation 

(1) |=/(^'-)' 

where the function /is expansible into a power-series in t and x, is included 
in most of the text books on differential equations. Methods have also been 
devised for proving the existence of an integral taking a given value ^ for 
t = T, when the function /satisfies much less stringent restrictions.* In many 
applications however the mere existence theorem is of little value unless ac- 
companied by a proof of the continuity and differentiability of the solution 
considered as a function of the initial constants t, ^. As the literature deal- 
ing with these latter questions is somewhat scattered and the proofs compli- 
cated, it is proposed in the following pages to collect and simplify as far as 
possible the results of the different writers. Solutions are defined reaching 
from boundary to boundary of the region R over which the properties of the 
function /are assumed provided It is closed, and the method of Peanot for 
proving the existence and continuity of the first partial derivatives with re- 
spect to the initial constants has been extended by a simple device to the 
higher derivatives. 

1. The existence of solutions. In order to make the existence 
proof suppose 1) that the function f{t, x) in equation (1) is continuous in 
the interior of a region! R of the ifa;-plane, and 2) that a positive constant 
k' can be found for any finite closed region R' interior to R, such that 

(2) \f{t,x)-f{t,x')\ S k'\x-x'\, 

* See for example, Osgood, Monatshefte fiir Mathematik, vol. 9, 1898. Only the continuity of 
f(t, X) is assumed. For otlier references, see Encyclopddie der mathematisehen Wissenschaften, 
11^ 4a, especially pp. 195, 200. 

t Atti di Torino, vol. 33, p. 9 ff. 

X In the following pages the word region means always a point set with interior points. A 
closed region Is one which contains all its limit points. 

(49) 
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where (t, x) and (t, x') are any two points for which the interval [x, x'] on the 
/-ordinate consists only of points in R'.* 

Consider the series of approximation functions :t 



(3) 






where (t, f) is a point interior to the region i? through which it is desired to 
pass a solution of equation (1). Any one x^"'> of these functions is surely 
well-defined on an interval 

(4) \t-T\^l 

for which the preceding function a;("-i> is continuous and the points (<, xC-^') 
all lie in the region R, and each function takes the initial value f when t = r. 
An interval (4) can be found by selecting a positive constant p so that all the 
points which satisfy the inequalities 

(5) /?i: \t-r\&p, \x-^\^p 

arc interior to the region R, and then taking I as the smaller of p and -—■ , 

Ml 

whore Mi is the maximum of the absolute value of fin the neighborhood (5) 

of the point (t, ^) . For then it follows from the relations 



.(«) 






^Mi\t-r\^p, (n=l, 2, ...), 



that all the points {t, a;'''), (<, a;'*)), etc., are in the neighborhood (5). 



194). 
then 



* This is the so-called Lipschitz condition (see EneyclopCtdie der Math. Wisa., IIA 4a, p. 

It will certainly be satisfied If / has a continuous partial derivative f^ = ^in S. For 

ox 



fit, X) —fit, x') = (x — x') fx{t, x + eix — z')), < e < 1 

provided that the interval [x, x'] of the «-ordinate is all in li', and the maximum of | fx\ in J?' 
■would be a constant k of the kind specified. 

t Compare Picard, Traite d'analyse, vol. 2, p. 301. 
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The series 

(6) a;=^ + (x('>-|)+ (xW-a;('))+ . . . 

is uniformly convergent* for all values of t on the interval (4). Namely 
from (2), 

(7) |a;W-a;("-i)| si f )/(<, a;(»-i)) -/(<,«(«-«)) jd/ 

I Jr 

g r^-ila;!"-!) - x(»-«) \dt , (n > 1), 

where ki is the constant in (2) belonging to the neighborhood ij!,- Let iVbe 
the maximum of \f(t, |) | on the interval (4) . Then 

and by successive application of the inequality (7), 

iV A;J|<-t|» 



|a;("> — xt"-!)! g 



^, 



The absolute value of any term of the series (6) is therefore less than the cor- 

JV 
responding term of the development of— jc*il'~^'— 1 j, which is aseries of 

Ki 

positive terms converging uniformly for ail values of t in the interval (4) . 
From this the uniform convergence of (6) follows. 

The continuous function x(t) defined by the series (6) on the interval (4), 
takes the value ^for t = t and satisfies the differential equation (1) . For since 
the region (5) is a closed one the function f is uniformly continuous, t i- e. 
for any positive e a positive S can be found such that 

\f(t', X') -f(t, x)\ <e, 

where (t',x') and (t,x) are any two points of B^ satisfying 

\t' - t\ < S, \x' -x\ < 8. 

P^urthermore on account of the uniform convergence of the series (6) a positive 
integer to can be found such that if n > m, 

|a;(»-») -x\ < B. 

* Here, as In what follows, the fundamental theorems concerning uniform convergence are 
assumed. See, e. g., Jordan, Conrs d'analyse, vol. 1, p. 310 ff. 

tA function which is coatiauous in a finite closed region, is uniformly continuous; see 
Jordan, loc. eit., p. 48. 
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Therefore also 



for h > m, and 



I f'\Ai,x^"-'^)-f(t,x)\dt \<d 
lira r/{t,x^«-^^)dC= ff{t,x)dt. 

n — cbJt Jt 



By letting m approach infinity on both sides of the equation (3) , it follows 
that 

(8) x = ^ + jf(t, x)dt, 

from which the statements made above are easily proved by putting t = t and 
by differentiation. 

A solution having been defined in the interval [t,t + i] , the same process 

can be applied again with the pointer + I, x(t + I) j as initial point, pi'ovided 

that this point is still within It, and the solution can be extended to an inter- 
val [t + ?, r + I + f] . If the region Jf is finite and closed the values of t 
which can be reached by such continuations must have an upper bound tj,&nd 

the points f<, x(t)\ approach a definite limit point as t approaches ^j. For 
if t' and i" both satisfy the incciuality 

where M is the maximum of |/| in the region Jl, then e<|uati(>n (8) shows that 

(9) \x{t') -x{l")\ 5 M \fJ -t"\ <€, 

and this is the necessary and sufficient condition that the function x{t) has a 
uni(]ue limiting value for < = /j.* The limit point fi,, a;(<i)j must be a bound- 
ary point of the region li. Otherwise it would be an interior point, and in 
that case /, would not be the upper bound of the <-va1ues which could be 
reached b^' contiimation. 



♦Suppose t', t", .... t'-\ . . . Imvc the limit «i and are less than h. By Jordan, loc. cit., p 
9, the correspomling values of x (lellue a iiniiiiie limit. Bat by (!>) this Is the limit which x ap- 
proiidics as t approaches I, over any set of valnes whatsoever. 
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On the other hand, if the region R is not finite and closed, it might 
happen that the solution could be extended over any ^-interval, or if t^ were 
finite perhaps the function x{t) would have no unique limit sX t = t^. In this 
case also, however, if <i is finite, and if x{t) approaches a unique limit as t ap- 
proaches <., then it can be proved as above that the pointful, a;(<j) ] must be a 

boundary point of R. 

The preceding discussion justifies then the following theorem : 

Through any poiiit (t,^) in the interior of the region R in which the 

function fis continuous and satisfies the Lipschitz condition (2), a solution of 

equation (1), denoted by 

(10) x = 4>{t,T,^), 

can be passed. The function <j) is continuous in t and defines points {t,x) in 
the interior of R for all values oft in an interval 

to < t < tj, 
in which the value r is included. 

If the region R is finite and closed, the interval [<o><i] is finite and can be 
chosen so that as t approaches t^ (or t^) , the points (t,x) of the solution ap- 
proach a unique limit point on the boundary of the region. 

If R is not finite and closed, the largest possible interval [t^, t{\ may be 
infinite, or if finite the solution may not have definite limit points at its extrem- 
ities <Q and ti. In case to (or <,) is finite aiid the solution approaches a unique 
limit point, then as before the limit point must belong to the boundary of R. 

For some applications it is important to notice that if P is a finite closed 
region entirely interior to R, a positive quantity X can be determined, such 
that the solution through any point (t, |) o/" P can be extended at least over the 
interval \t — r\ S X. This follows because a constant p can be chosen inde- 
pendent of the position of (t, |) in P and so that the corresponding region (5) 
is entirely interior to R. If M is the maximum of |/| in the region consist- 
ing of all the regions (5) about points of P, then according to the preceding 

discussion the smaller of p and ^ will be a constant X with the desired pro- 
perty. 

2. Uniqueness of the solution.* Considera solution 6" (x = a;(<)'\ 
of equation (1), defined and interior to the region R for all values of t in an 
•Compare the methods nsed in §§2, 'A, i with Peano, loc cit. 
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interval [f, ^']. A positive constant p can be so selected that all points of the 
neighborhood S^ of 8 indicated by the dotted lines in the figure, are also in- 
terior to a. Let sly = i/(t)jhe another solution such that the value t = r. 




and therefore all values in a certain interval T including t, define points (t, y) 
interior to the region R. Then the difference (x — y) satisfies the equation 

and the forward derivative* of |a; — yj satisfies the inequality 



(11) 



\d\x — y\ 



dt 






dt 



on account of the condition (2) , where k^ is the value of k belonging to the 
region R^. Consequently 

d\x - y\e-^<f—^ ^ ^ d|a;-y|e*'(«-^> ^ ^ 
dt - "' dt ^ "• 

By using the former inequality for / > t and the latter for < < t, it follows that 

(12) |»_y| ^ If _,|e*.l<-»l 

for all values of t on the interval T, where | and t] are the values of x and y 

*Th« forward derlratire is used becaose at points wliere a function u(() vanisties, tlie de- 
rivative of I tt I s= + Ju* is discontinuons. At such values of t, however, the forward deriva- 
tive exists, and by the mean-value theorem Is 

^' _ lim 1 
dt ~ »=« h 

where A is supposed to approach zero from the positive side. At other pointstbeflrst partof (11) 
is easily proved. 



-WW. 
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for t = r.* It follows, then, that if 8 and S intersect at any point t = t, ^ = r), 
the value of \x — y\ is less than or equal to zero, and x must be identically 
equal to y for all values of t in the common interval adjoining t in which the 
two solutions are defined. 

Only one solution of equation (1) can pass through a given point (t, |) 
interior to H. The set (10) includes, therefore, all the solutions in the region. 
3. Continuity of the solutions with respect to the initial 
values. With the help of the inequality (12) it can be shown that the solu- 
tion (10) through the point (t, f) can be surrounded by solutions correspond- 
ing to neighboring initial values (t + At, | + Af), and that these go over 
continuously into the former when At and Af approach zero. Let S denote 
the solution (10), and suppose that an interval [<', <"] and a neighborhood i?2 
of the kind described in §2 have been found. At and Af are to be so chosen 
that the point (t + At, f + Af ) lies in B^, and S will be used to denote the 
solution 

X = <f>(t, T + At, I + A|) 

through that point. Then from equation (8) 

T-f At 

|^(t + At, t+ At, | + A|)-<^(t + At, t, 1)1= A?-//(<,0)c 

^\A^\+M,\Ar\, 
where M^ is the maximum of | f\ in the neighborhood i?2' This gives the dif- 
ference of the values of x for S and 8 a.t t — r + At, so that when r is replaced 
by T + At, the inequality (12) becomes 

(13) |<^(<,t + At, 1+ A|)-<^(<,T,f)|^{|Af|+if2|AT||e*^i'— M. 

Suppose now that At and Af are still further restricted so that 

(14) j |A|| + Mj|At| je*^i«"-«'l < p. 

Then the solution 8 through (t + At, f + Af ) can be continued over the whole 

*Ifa function f(t) is continuous and has a forward derivative f(t) ^0 on the interval from 
tt, to ti (fo < «i), thenf(t(,) S/(«i). For, the function 

where /S is any negative constant, has a negative forward derivative. If ^(io) were less than 
(t>(,ti), the function (p would attain its minimum at some value t' less than ti. At that point the 
forward derivative would necessarily be S 0, contrary to the hypothesis. Consequently 

*(«o) - *(«i) =/(«o) —f(ti) — /3(«, - (o) > 
for every negative p, which proves the theorem stated for/(«). 
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interval t'f. For by the results of § 1 , the process of continuation can be car- 
ried on until the boundary of i?2 is reached, and on account of the inequalities 
(13) and (14) the only points attained by such continuations between the or- 
dinates t = t' and t = f are in the interior of R^. The ordinates < = i' and 
t = I" are therefore the only boundaries which the solution 8 can cross. 

From these results it follows that '\i t, t -\- -At are any two values in the 
interval \_t', t"], then by applying (8) and (13), 

\A4>\^\<t>{t + At, t+At, ^ + A^) - 4,(t, t, ^) \ 

^\<f>(t + At, t+At, f + A?) -</)(«, t + At, 1 + A^) | 

+ \4>{t,T + AT, ^ + A^)-<t>{t,T,^)\ 

^Mi\At\+ \\A^\+ M^IAtI |e*!(«-^-^^). 

Therefore the function (10) is continuous for all values of t, t, ^ defining 
points of the solution interior to the region R. For, as has been indicated, an 
interval [f, <"] including the values t and t can always be selected, and 
At, A^ restricted, so that the solutions S are defined on the interval and sat- 
isfy the last inequality. 

4. Partial derivatives with respect to the constants of inte- 

gration. In order to prove the existence of the partial derivatives ^ , -^ 
of the function (10) at a point t, t, f, the additional assumption is made that 
the derivative /a; = r=^ exists and is continuous in the interior of the region R. 

For the two solutions 8, 6* of §3, it follows from equation (1) that 

(15) ^ = »A*, 

where 

A4> = <l>(t, t + At, I + A|) - <f>(t, T, ^) , 

^ f(t,<l> + A4>)-f(t,<f>) 

""- A^ 

In these expressions t and ^ are considered fixed, for the moment, while t, At, 
A| are variable. The coefficient a is a continuous function of t. At, A^ when 
t is on the interval [f, <"] , and At, A| satisfy (14) .* This follows from §3 for 
any values of the arguments for which A<f> i^ 0, and with the help of the expres- 
sion 

a^f^(t,^ + eAi,), 0<^<1, 

♦See a note by Hadamard, H'ull. de la societe math., vol. 28 (1900), p. 64. 
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for any values for which A<^ = 0, on account of the continuity of the deriva- 
tive fx- The general integral of equation (15) is therefore 

A0 = ce 

When At = 0, the constant c has the value Af, and it follows easily that 

AA 
the quotient -^ has the limit 

(16) I ^ //.<'.*)'« ^. 

which is a continuous function of t, r, ^ on account of the continuity of <f> 
and^. Similarly when Af = 0, 

c = A<^|'=^ = r f{t, <f> + A(f>)dt 

Jr+ At 

= - At/(t + Mt, 4>(r + ^At, t + At, |)^ , < ^ < 1, 

by (8) and the mean value theorem for a definite integral. At the limit there- 
fore 

(H) f^ = -/(.,l)/'^^^''"" 

Differentiation of equation (8) twice for t shows that it the function /has 

df . . 3»^ 

the derivative ^ continuous in the interior of B, then the derivative -^ ex- 
ists and is also continuous. If furthermore all the second derivatives of/ ex- 
cept perhaps ^^ are continuous, the equations (16) and (17) can be differen- 

tiated, and ^ will have all the derivatives of the second order. A simple 
induction leads to the theorem : 

When the function f(t, x) has all. t<.s derivatives of (he (n — 1)" order, 

9"/' 
and all of the n'" order except perhaps -^ , continuous in the interior of the re- 
gion R, then the solution (10) has continuous derivatives of the n"' order with 
reajject to t, t, ^ at any set of values (t, t, f) defining a point {t, x) of the solu- 
tion interior to M. 



* See Beudixon, Jinll. de la socii'U math., vol. 24 (1896), p. 220. 
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5. Complex numbers.* The results of the preceding sections can 
be readily extended to a system of differential equations by means of complex 
numbers and matrices. The necessary definitions follow : 

1) x is a symbol for the set (xj, asj, • • • , x,.) ; 

2) mod X = y/xi + x^ + . . . + xl'. ; 

3) cfgyifx, §y, (i=l,2, . . .r); 

4) X ± y are the complex numbers (Xj ± yi, x.^ ± y.j, • • • , x,. ± y^) ; 

r 

5) xy is the sum 2 x^y^ ; 

« = i 

6) kx is the complex number (kxi, kx^, • • • , kx^) if jfc is a simple 
constant ; 

/dx 
xdt and -=- are the 

complex numbers ( jxydt, fx^dt, ■ ■ . , fxrdt\ and ("^ , ^ , . . . , ^\ 

respectively ; 

8) |x| is the complex number (|xi|, \x^\, • • • , |Xr|) ; 

9) A denotes a matrix of elements O/j (*>/= 1> 2, •••,»•) ; 

10) Ax is the complex number of which the elements are a,iXi + ai^.j. 
+ • ■ • + a,;^^ (iz:^l,2, ■ . ■ ,r); 

11) A = Bi{ a,j = h^j (i, J =.1,2, ... ,r); 

12) AB is the matrix C whose elements are Cy = ajfiyj + a,-.26.^- 
+ • • • + ttirbrj (i,j = 1,2,...,r); 

13) —- is the matrix whose elements are -j^ (h J = 1, 2, . . . r) ; 

14) -^, whenx and ^ are two complex numbers such that the ele- 
ments of X are functions of the elements of ^, is the matrix of derivatives 
|| (i,J^l,2,...r). 

Besides these definitions the following three theorems are needed : 
a) mod(x ± y) S mod x + mody ; 



♦For the application of complex numbers to linear equations, see Pcano, Math. Annalen, vol. 
32 (1888), p. 450. In the article by Peano cited previously, they are used without explanation 
for any system of equations in the so-called normal form. 
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b) if the elements of x are functions of t having continuous deriva- 
tives, then mod x has a forward derivative and 



I (2 modx 



c) mod 



/'•" ^ 1/ 



dt 
modxdtl 



& mod 



dx 
dt 



To prove a) it is only necessary to apply the identity 

(18) i^' iy\=(ix^t\\ i (xnjj - xj7jty, 

after squaring both sides of the inequality. To prove b) two kinds of ^-points 
must be considered, according as x ?i or x = 0. For the former [sec 5)] 






rfmodx 


= 


dx 
''di 
modx 


dt 



on account of (18) when y 
the mean value theorem is 



dx 
dt 



; and for the latter the forward derivative by 



^ = is I V,!,[*"^('-^")J' = - V,!,(w) • 

(see footnote p. 54) where k approaches zero from the positive side. From 
6) it follows further that 



Therefore 



-r mod / xdt ^ modx = t- / i 



mod xdt. 



1 
dt 



I mod xdt — I mod xdt I ^ 0, -r ] mod xdt + I mod xdt I S 0, 



and by using the former for t > t and the latter for < < t, it is seen that <:) 
must be true on account of the footnote on p. 55. 

6. Systems of differential equations. By means of the definitions 
of §5 a system of differential equations 



(19) 



dXi _ „. . 

—fl —JiVt ^l> X^t • • •> Xr), 



= 1,2, 



,'•). 
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can be represented by a single equation of tlie form (1), where however x and 
fare complex with the elements Xi and^ respectively. It is presupposed that 
the functions fi are continuous in the interior of a region i? of r + 1 dimen- 
sions, and that a positive constant V exists for any finite closed region B! in- 
terior to i?, such that 

(20) mod[/(^ X) -fit, y)] g U mod (a; - y) 

whenever (<, cc), (<, y), and all the points (t, z) for which the elements Zi are 
in the interval [a;,-, y,], lie in the region B! .* 

If ^ is a complex number and (t, f ) is any point interior to R, a positive 
constant jO can always be found such that all points (t, x) satisfying 

i?i : \t — t\ S: p, mod (x — ^) ^ p, 

are also interior points. The set of approximation functions for the system 
(19) is defined by equations of the same form as equations (3) but with 
^(>t-i)^ xW, ^, and/complex. li t remains on the interval 

(21) \t-T\ g?,t 

where I is the smaller of p and -^ ,Mi being the maximum of mod ft in the 

Ml 

*The analogue to the Lipschitz condition of §1 would be, in the notation of complex num- 
bers [see §4, 8), 10), and 3)], 

(21) \At,x)-f(t,y)\ ^ K \x-y\, 

■vvhere K is a matrix of positive constants kij. The condition given In the text is a consequence 
of this, as can be showu bj' taking the moduli of each side of the equation (21) and applying the 
relation 



»,.y=i 






which holds for any quadratic form because j -t'l j ^j S «j 



'' Ni.' 



^1. 



In case the elements fi have continuous partial derivatives with respect to the variables x 

the value of i„ can be taken equal to the maximum of the absolute value of /■», = — (i, 1 = 1, 

■' d%j ^ 
2, . . . , r) in the region B. For 

/•,(«. J-i, . .. ,*,.) —/■(«, 2/1 2/,) = .2 (xj — yj)/tj(^t,yi + e^ixi— yi),. . . ,x, + 0,(x^ — Vr)), 

where the d's are all between zero and unity. 

fit should be noticed that if the system of equations is linear in the xi's this restriction on 
t is unnecessary. The approximation functions are well-deflned on any interval for which the 
coefficients of the Xi's are continuous functions of f. (See §6.) 

Jin many cases the proofs for a. system of equations can be worded as for a single equation- 
provided that absolute values are replaced by moduli as defined in §4. 
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region lii, then by means of §5 c) it follows step by step that the functions a;'"' 
satisfy the relations 

mod(a;('0 - ^) ^ Mi\t - t\ ^ p, 

and are well-defined and continuous in the interval (21). 

In a manner analogous to that of § 1 , it can be shown with the help of (20) 
and §5 c) that the terms of the series 

(22) mod(a!(" - t) + mod(a;W - «<») + mod(«W - «'«) + • • • 

are less than the corresponding terms of 

^|^iK-T| + 2j +— 3j +...J, 

where JV is the maximum of mod f(t, f) in the interval (21), and the series 

(22) converges therefore uniformly on the interval (21) . The r series defined 

by 

(23) x = ^ + (x(i) - f) + (xW - xW) + . . . 
are all convergent, and as in § 1 , 

(24) x = ^+ly{t,x)dt. 

Hence the series (23) define r functions Xi(t) which are continuous on the in- 
terval \t — t\ ^ I, have the values ^, for t = r, and satisfy the differential 
equations (19). 

These solutions can be continued to the right and left of the interval 
[t — 7, T + if] , and the values of t which can be reached by such extensions 
have a lower bound t^ and an upper bound ty as before, provided the region 
R is finite and closed. 

Through any point (t, fj, ^j, •••,!,) interior to the region R in which 
the properties of the functions f are dejined, there passes a. .fystfo of solutionis 
of the equations (19), denoted by 

(25) Xi = 4,i(t, T, ^1, ^2, . . • , ^r) {1=1,2, ... ,r), 
or, in the notation of complex iiumbers, 

x = <f>(t, T, I). 

The functions <f>i are continuous in t and define points (t, x^, x^, . • . , x„) in 
the interior of R for all values oft in an interval 

<o < < < <i 
including the value t. 
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The statements made in §1 concerning the behavior of the solution in the 
vicinity of t^ and <i apply also to the solutions (25). 

7. Uniqueness and continuity for a system of equations. 
The region H^ (compare §2) in the neighborhood of a system S of solutions 
X = x{t) which are defined and are interior to i? in an interval [<', <"] , consists 
of all the points (t, x) satisfying the inequalities 

t' ^t ^ t", mod[a; - «(<)] ^ p. 

For a properly chosen p, R^ will be entirely interior to i?. K S(y = y{t)^ 

is another system of solutions interior to H.^ for t = t, then by §5 h) and (20), 

|mod(a; - y) S mod^i^^^ = mod[/(^ x) -/( t, y)] 
^ A;2mod(a; — y), 
where Xvj belongs to li^ . As in §2, it follows that 

(26) mod(a; - y) ^ e*-=i'-'i • mod (| - t]), 

where ^ and tj are complex numbers representing the values of a; and y at 
t = r. But mod (x — y) is necessarily S 0. Hence on/y owe system of solu- 
tions of equations (19) can pass through a given point (t, |j, I.^, • • • , |„) 
interior to R. The functions (2o) represent therefore all systems of solutions 
in the region. 

Let the system x — x(t) be replaced by the system of solutions (25) , and 
let Ji^ be so chosen that the interval [<', <"] includes the value t. If At and a 
complex number A^ = (A^i, Afj, • • • , A|„) are such that the point (t + At, 
I + A^) is interior to i?^, then from (26) when t is replaced by t + At, 

mod[<^(<,T + AT,^+A^) -<^(/,t, I)] 

g ei-,l«-T-Ar|,„od [1 + AJ-^(t + At, t, I)] ; 
or by means of (24), 

mod[<^(<, T + A t, I + A|) - <f>{t, T, ^)] 

Se'-^i'-'-^'lmodfA^- T ''/(<, ^)rf<l 

But by §4 a) and c), this becomes 

(27) mod[</.(<, T + At, I + AO - <^(<, r, f)] 

S e'>,l'-'--^" [mod A| + J/jIAtI], 
where M.^ is the maximum of mod/in i?.,- The inequality (27) holds through- 
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out the interval T including t, in which the solution 8, 

X = ^{t, T + At, ^ + Af) , 

through the point (t + At, f + A^) is defined and interior to R^. 

Choose now At and A^ so that 
(28) e*,i«"-''i [modA|+ ilfj I At I ] < p. 

The solution /8'can then be continued over the whole interval \t', <"] because, 
on account of (28) , S remains interior to i?2 for all values of t between t' and t" . 
Furthermore if t and t + A< are both in [<', <"] , 

mod A</) = mod [<^(< + A<, t + At, | + A^) — ^{t, t, f)] 

^Mi\M\ + e*=l'-'-^H[modA| + il/j|AT|], 

by means of §5 a), (24), and (27), in a manner similar to that of §3. It is 
readily seen that the solutions (25) are continuous for all values of I and 
the initial values (t, ^i, f.j, • • • , fr) which define points of the solution in- 
terior to the region R. 

8. Systems of linear equations. The system of linear equations 

(^9) II =i^^^^ (i=l,2,...,r), 

can be represented in the notation of complex numbers by the single equation 

dz 

dt='^' 
[see §5, 7), 9), 10)], where z is complex with elements Zf, and A is the 
matrix of elements a,-, (i,y = 1, 2, • • • , r). The coefficients ay are supposed 
to be continuous functions of t on an interval <o ^ < = h- By the methods ot 
§§6, 7* a matrix E can be determined, of xvhich the columns are r systems of 
solutions of the equations (29), taking respectively/ the columns of the matrix 

1, 0, 0, . . . , 

0, 1, 0, . . . , 

(30) «' <>• 1. • • • , 

0, 0, (),..., 1 



*The region R consists of all points («, z) for which („ S « g «, and z is arbitrary. The 
constant* can be taken as r times the maximum of the absolute values of the a'J in the interval 
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as initial values for t = t. The elements of E are continuous functions of t 
and T lohen these values lie in the interval [4, <i] • 

From §5, 12) , and the fact that the columns of E are systems of sohitions 
of equations (29), it follows that 

Consider the set of functions 

2 = -£■?, 
where 5^ is a complex constant. By differentiation and application of (31) 
[see §5, 7), 13), 10), 12)], 

dz dE ^ . _^ . 

Hence the system of integrals of equations (29) talcing any given initial values 
^for t = r, and which by ^1 is unique, can be expressed as the product of the 
matrix E and the initial valves ^, 

z^E^. 

In the next section it will be important to consider the solutions of linear 
equations in which the coefficients a^j involve a number of parameters besides 
t. Use will be made of the following 

Auxiliary Theorem : Jf the coefficients ay of the equations (29) involve 

also s parameters a — (aj, a.^, ■ ■ ■ , a^), and are continuous functions of them 

in a finite closed region P, then the elements of the matrix E are continuous 

functions oft, t, a for all values of t, t, a such that t and t are in the interval 

[<o, <i] and a in P. 

The proof is easily made by examining the set of approximation functions 
of the form (3). Each of them is defined on the whole interval [<o, 'i] (see 
footnote p. 60), and is a continuous function of t, r, a. If the constant ^of 
§6 is taken as the maximum of mod/(<, f) = mod A^ for all values of t in 
[to, <i] and of a in P, then the convergence proof is independent of the param- 
eters a, and the series corresponding to (23) are uniformly convergent 
series of continuous functions of t, t, a; which proves the theorem.* 

9. Partial derivatives of a system of solutions. The existence 
of the first partial derivatives of the solutions (25) with respect to the initial 

*See .Jordan, loc. eit. p. 311. 
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constants can be proved, provided that the functions j^(<, Xj, x^, • • • , x„) 
have continuous first partial derivatives with respect to the a;'s in the interior 
of the region i?. Let t, r, | define a point of (25) interior to R, and suppose 
74 constructed as in §7 so that the interval [f, <"] includes the values t and t. 
The difference of the two solutions <f), (f> + A<f> corresponding to (t, |) and 
(t + At, I + A^) respectively, satisfies the equation 

(32) ^ =/(<, <j> + Acf>) -f(t, <!>) = AA<f>, 
where the elements of the matrix A are determined by the equations 

— /•(^ <t>l, ■ • ■ ' 'f>j-U <t>jy 'fc + l + ^ff^j + l, • • • , <f>r + ^<t>r)- 

If T and I are considered fixed for the moment, it can be shown (compare with 
§4) that the coeiBcients a/j are continuous functions of t, At, A^, provided that 

(33) t'^t^t", |At| = |t-t| ^ S, modA^ = mod(^ - I) ^ S, 

wliere S is so chosen that the point (t + At, | + A^) is in H.^ and the inequality 
(28) is satisfied. This follows from the fact that the solutions (f> + A<f> corre- 
sponding to such values of At, A^ arc defined and interior toi?2 in the interval 
[(', t"], and according to §7 the functions <f> are continuous. 

By the auxiliary theorem of §8 a matrix U exists for the equations (32), 
with elements which take the initial values (30) at < = t, and which are con- 
tinuous functions of t, At, Af in the region (33) . Suppose that A^j is diff"erent 

from zero, while At= Afj = Afj =...=: Af^ = 0. The r quotients -r-^ [see 

§5, (!)] are a system of integrals of equations (32) with the initial values 
(1 , 0, 0, • • • , 0) for t = T, and therefore coincide with the first column of U. 

A similar reasoning applies to the other systems -^ (j=2, ■ ■ ■ , r) for 

which At and all the elements of A^ except A^^ are set equal to zero. Since 
the elements of U are continuous functions, it follows for each system that 

the liniits of the ([uotionts -—■ as A^^- approaches zero exist, and at the limit 
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(34) ?| = ^|A^ = Af=«, 

where ;^ is the matrix of derivatives ^ (i,j = 1, 2, ■ ■ ■ , r)- 

So far only the existence of tlie derivatives (34) has been proved, but it 
can be readily shown that they are also continuous at the given point t, t, |- 
By substituting a; = <^ in the ecjuations (19) and deriving successively for 

^, (i = 1, 2, ■ • ■ , »•), it follows that the columns of the matrix 5? are solu- 
tions of the equations 

dz , 

(35) ^ =/.., 

where/,, denotes the matrix of elements -^ (i, J = 1, 2, ■ • • , r) in which 

OXj 

the arguments are t, <f>(t, t, |) . Equations (32) go over into these equations 
when At = A^ =r 0. If the rand | are replaced by t = t + At and f = f + A|, 

the cootlicients ~- will be continuous functions of t, t, | in the region (33). 

The elements of the matrix of solutions -^ with the initial values (30) at < = t, 

arc therefore continuous functions of ^, t, | in the region (33), in particular at 
the point (t, T, I). 

When At jt 0, A| = 0, the value of A(f> at < = t is found from equa- 
tion (24) to be 



(3()) A.^ 



t = T 



•'t + At 



+ A(l>)dt= - At . /, 



whore the elements of the complex number f are 



./;. (^T + eAr, 4>(t + OAr, T + At, f) V < 0, < 1, (i =1, 2, ■ ■ ■ , r) . 

Accordingly the ((uotients -j^ , which form a system of solutions of equatioi 
(32) with initial values at< = t defined by (3()), have the values (see §8) 

A7 = - ^/- 
At the limit for At — 0, 
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where the right member is the product of the matrix -^ and the complex num- 

ber/(T, f) [see §5, 6) and 10)]. 

It has been proved therefore that wAen the functions fi(t, x^, x^, ■ ■ ■ , x„) 
have continuous first partial derivatives with respect to the X{S in the interior 
of the region li, then the solutions (25) have first partial derivatives which are 
continuous for all values t, t, f defining points of the solution interior to li. 

The existence of higher derivatives at a point t, t, | when further as- 
sumptions are made upon the functions f, can be proved by means of the 
2r + 1 equations, 

where the arguments of f. are t, (f)(t t, |). The right members are functions 
of t, T, I, z, and arc continuous when t, r, | satisfy the conditions (;5;$) wliat- 
evcr the vahic of the arguments z. If the f have continuous second partial 
derivatives with respect to the aj/s, then according to the last theorem the right 
members of (38) have continuous first partial derivatives for t, |, z, and their 
solutions 

are therefore difTerentiable for the initial values t, f. If furthermore the f 
have continuous first partial derivatives with respect to t, then on account of 

the ditlercntiability of the functions ^ and/-, the expressions (37) can be 

differentiated once, and <t> twice, for t. From (24) the functions <p can also 
be dillercntiated twice for t. In general the following theorem can be proved : 

If the functions f(t, Xi, x.^, • • • , x,.) have all the partial derivatives of 
order n — \, and all of order n except perhaps the n"' derivatives with respect 
to t alone, continuous in the region J{,then the sohtiions (25) of equations (19) 
liave continuous partial derivatives of the 9i"' order for all values of I, t, ^ 
defining points of the solution (25) interior to It. 

For suppose that it holds for n, and that the functions f have all deriva- 
tives of the ?i'" order and all of order n + 1 except perhaps the (?i + 1)" de- 
rivatives with respect to t. Then according to the hypothesis the functions 
have continuous «"' derivatives, and, on account of the existence of the 
(n + 1)" derivatives of the functions y), the right members of equations (3«) 
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can be differentiated »t times for <, t, |, and z. The solutions ~ and the ex- 

pressions (i??) can therefore also be differcjitiatcd n times. Since the func- 
tions fi have continuous )i"' derivatives with respect to t alone, equation (24) 
shows that the functions ^,- can be differentiated n + 1 times for t alone. The 
induction is therefore complete, and the theorem holds for any value of n. 

The Univeusity of Missouhi, 
Columbia, Mo. 



